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ABSTRACT 

The anomalous chemical abundances and the structure of the Edgewood-Kuiper belt observed in the solar 
system constrain the initial mass and radius of the star cluster in which the sun was born to M ~ 500 to 3000 Mq 
and R ~ 1 to 3pc. When the cluster dissolved the siblings of the sun dispersed through the galaxy, but they 
remained on a similar orbit around the Galactic center. Today these stars hide among the field stars, but 10 to 60 
of them are still present within a distance of ^ 100 pc. These siblings of the sun can be identified by accurate 
measurements of their chemical abundances, positions and their velocities. Finding even a few will strongly 
constrain the parameters of the parental star cluster and the location in the Galaxy where we were born. 



1. INTRODUCTION 

It is c ommonly a ccepted that stars like the sun are born in 
clusters (iLada et al.L [1993). The cluster in which the sun was 
born is long gone and the sun's siblings are by now spread 
over the Galaxy. The structure of the hot Edgeworth-Kuiper 
belt objects provide evidence of this, as this population can 
be r eproduced by a relatively n earby encounter with another 
star jMorbidelli & Levisonl 12004) . Such an encounter is ex- 
pected to have oc curred in the early history of the Solar system 
(lMalhotraLl2008h . The existence of a well organized planetary 
system, however, indicates that the parental cluster cannot have 
been very dense as otherwise a nearby passing star would also 
have perturbed the orbits of the planets. 

Additional evidence for the sun's dynamic history comes 
from the discovery of radioactive isotopes and their decay prod- 
ucts in the proto-solar nebula (Hester et al., 2004), which is ex- 
plained by a superno va explosion within L6 pc of the infant sun 
jLoonev et al.L l2006l) . The combination of arguments enables 
us to estimate the mass and size of the star cluster in which the 
sun was born. We follow the orbital evolution of these stars 
through the Galaxy for the lifetime of the sun and conclude that 
at least 1% (10 to 60) of the sun's siblings should still be present 
within lOOpc, and more than 10% should be within about a 
kpc along the orbit of the solar system in the Ga lactic poten- 
tial. With the Gaia satellite dPerrvman et all 1200 1[) and ground- 
based searches the radial velocity and distance to the majority 
of the lost siblings will be determined, and the proper motion 
will be measured. These constrain the orbit of the proto-solar 
cluster and enable us to accurately determine the evolution of 
the Galactic potential and the birth place of the sun. 

2. THE PARENTAL STAR CLUSTER 

The sun and its eight planets were born about 4.57 billion 
years ago ( Bonanno et a 11 12002 ). probably in a star cluster, and 
it is in due time that our location in the Galaxy became so deso- 
late. Evidence for this dynamic past comes from meteorite fos- 
sil records, where the presence of short-lived radioactive iso- 
topes in primitive meteorites indicate that the 1.8 Myr young 
sun was polluted by a supernova explosion of a star about 15 to 
25 times m ore massive than the sun within a distance of 0.02 to 
1.6 parsec ( Loonev et al.L l2006l) . Such a massive star lives for 
6 to 12 Myr before it sheds the majority of its mass in a super- 



nova explosion. This massive star must have formed some 4.2 
to 10.2 Myr earlier than the sun. Such range in the distribution 
of stellar ages is also observed in the Orion Nebula where mas- 
sive stars also tend to be a few to ~ 10 Myr younger than the 
low-mass stars (Palla & Stabler, 1999). 

The presence of a massive star close to the infant sun 
puts interesting constraints on its birth environment. Today, 
such a massive star is not even present within a distance of 
100 parsec. By a dopting a standard initial mass function 
dKroupa & Weidnen , .2003) about 1 in 400 stars is sufficiently 
massive to experience a supernova. Stars of to = 15 to 
25 M0 are less common and would require a star cluster of 
M > 500 M© (IWeidner & Kroupa, 2004). 

Massive stars in a cluster tend to sink to the center in a frac- 
tion oc 1/m of the two-body relaxation time scale (irix)- If 
trix ^ 300 Myr even the most massive stars are unlikely to have 
reached the cluster center by the time they explode. Massive 
stars in clusters with smaller tiix tend to populate the central 
region by the time they explode in supernovae, polluting the 
nearby young stars and proto planetary disks in the process 
(IWielen et all Il996h . Low mass stars, like the sun, are not 
strongly affected by mass segregation, and if such a star hap- 
pens to be in the cluster core at the moment one star explodes it 
is likely to be still around when a subsequent supernova occurs. 

We quantify this by performing a number of iV-body simula- 
tions in which we initialized star clusters wit h a mass function 
and a Plummer sphere density distribution (iPlumme 
with radius R in parsec to track the number of supernovae that 
occurred within 1.6pc of sun-like main-sequence stars. It turns 
out that the number of supernovae with one and the same star 
of 0.8 to 1.2 Mq scales as iVncarbySN ~ 0M3M/R. 

A star cluster is a dynamic environment in which fre quent en- 
count ers effectively ionize planetary systems (Malmberg et al.L 
I2007h . The passage of a star within a distance of renc = 100 
Astronomical Units (AU) can easily destroy the young proto- 
planetary nebula, leaving the sun without any material to form a 
regular planetary system (lHillslll985llKobavashi & Idall2001h . 
If the planetary system is already formed such an encounter 
may induce high eccentricities and inclinations of their or- 
bits, again resulting in strong perturbations which event ually 
leaves the sun with a reduced planetary system (Spurzem et al.L 
I2OO6I) . The high eccentricities and incUnations observed in sev- 
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eral extra-solar planetary systems may be produced by such a 
relatively nearby encounter. The outer region of the solar sys- 
tem shows evidenc e for a relatively nearby encounter but not 
as close as 100 AU jlda et al.. 2000), an encounter at a distance 
within Tone — 10'^ AU would be sufficient to explain the ob- 
serv ed hot Kuiper belt objects, which have highl y inclined or- 
bits (iMorris & 0'NeillLll988tlHester et al.Ll2004 . 

The requirement that the sun has not encountered another star 
at Tone ^100 AU, but is likely to have experienced an encounter 
within some r^nc ^ 10'^ AU constrains the initial cluster den- 
sity. We further require that the cluster survives long enough to 
warrant this encounter to occur before it dissolves. The dissolu- 
tion time of a star cluster in the MiUcy W ay Galaxy follows the 
empirical relation jLamers et al.L l2005h tdiss ^ 2.3Myr A/° '^. 
These constraints can be expressed in terms of cluster radius R 
(in parsec) and mass M (in Mq), which we present in fig.[T] 
With estimates for the mass and the size of the parental star 
cluster we calculate that the velocity dispersion of the clus- 
ter must have been of the order of a km/s, which is compa- 
rable to the velocity dispersion in known young star clusters 
(iLada&LadaLlIool . 

3. THE BIRTHPLACE OF THE SUN 

The orbit of the solar system in the Galactic potential is al- 
most circular, more or less in the disk. The orbital velocity is 
considerably higher than the velocity dispersion of the stars in 
the parental cluster The cluster members are therefore unlikely 
to have drifted very far from the orbit of the sun around the 
Galactic center, but they may be at completely different loca- 
tions along this orbit. The preservation of phase space in the 
dissolution of star clusters is used to study the merger history 
of the Galaxy jHelmi et al.L [1999 ). but in this case we use the 
argument to find the lost siblings of the sun. 

At the moment the sun is located at a distance of about 
8.5 kpc from the Galactic center in the plane of the disk. The 
velocity is 10.1 km/s directed towards the Galactic center and 
7.5 km/s away from the plane. The velocity perpendicular to 
the line-of-sight of the Galactic center is vt = 235.5 km/s in the 
rotation direction of the Galaxy, which is slightly (^ 15km/s, 
since c± 220km/s) hig her than the local orbital velocity 
jKarachentsev & Makaro\II 19961) . 

By calculating the orbit of the sun backwards in time, which 
can be achieved by inverting the current velocity vector, and in- 
tegrating the equations of motion in the potential of the Galaxy 
for 4.6 Gyr, we can calculate the location in the Galaxy where 
the sun was born. With the parameters giv en and using a model 
for the Galactic potential ( Paczvnskilll990 ) it turns out that the 
sun has been orbiting the Galactic center some 27 times since its 
formation, and was born in the (5-quadrant at -1.39kpc, 9.34kpc 
and 25.3pc along the galactic x, y and z-axis and with a veloc- 
ity vector of -207, -48.2 and -6.72km/s (see fig.©. This birth- 
place is ~ 2.8kpc further out than the estimated distance to the 
Gala ctic center based on the relatively hi gh metalicity of the 
sun jWielen et al.L 119961; IWielen & WUsonl 119971) . but part of 
this metalicity argument could be explained by the early pol- 
lution of a supernova near the proto-solar nebula. Part (up to 
about half) of the other cluster members are also expected to be 
polluted by the same supernova, and their anom alous chemical 
chara cteristics can be used for the identification dTolstov et al.L 
120041) . 

It seems likely that the orbits of the sun and it siblings have 
been deflected by some small-angle scatterings on their long 
journey through Galactic disk. The radial distance over which 



the sun may ha ve been migrated by suc h scatterings can be as 
large as 2 kpc dFuchs & Wielenl Il987l) . which interestingy is 
comparable to the range in radial distance for the orbit of the 
sun, as we present in Fig.|4] The radial velocity induced by 
such scatterings affect the orbit of all stars in the parent clus- 
ter and the effect on the number of nearby stars today is rather 
small. We quantify this by varying the current tangent veloc- 
ity of the sun vt and recompute the orbits of its siblings in 
the potential of the Galaxy assumig that the cluster dissolved 
4.3 Gyr ago. The fraction of stars of the parent cluster that re- 
main within lOOpc can then be expressed as a linear relation 
/(r < lOOpc) ~ 0.032 - 0.014('i;t/wc), and the fraction of 
stars within Ikpc becomes f{r < Ikpc) ~ 0.39 — 0.16(ut/wc)- 
These relations break down then vt is either so small that the 
stars enter the Galactic bulge on their orbit, in which case the 
scattering becomes severe, or if vt becomes sufficiently high 
that the stars escape the Galaxy. 

The uncertainty in the current coordinates of the sun in the 
Galaxy will be magnified by computing its orbit back with 
time, making the location of birth uncertain. It turns out, how- 
ever, that the exact place and velocity of birth are not crucial 
for our discussion, as the stars that were born together with 
the sun will be around wherever the sun was born. It is the 
orbits of those stars for which we integrate the equations of 
motion in time through the Galaxy to understand the disper- 
sion in their position and velocity at a later time. The adopted 
model for the Galactic potential does not include spiral arms 
or local perturbations like star clusters and molecular clouds, 
which can have an appreciable effect on the calculated t rajec- 
tories dFuchs & Wielenlll987l:l0uillen & Minchevil2005h . The 
uncertainties introduced by the model, however, have a minor 
effect on the relative position and velocities as all stars remain 
relatively close together in phase space throughout the integra- 
tion and are therefore affected by variations in the potential in 
a similar fashion. 

4. FINDING THE LOST SIBLINGS OF THE SUN 

We persue by constructing model star clusters that mimic the 
one in which the sun was born. These star clusters are assumed 
to dissolve along the trajectory of the sun in the Galaxy, shed- 
ding their stars in radial orbits with the escape speed of the clus- 
ter. In that way we are able to calculate the fraction of the stars 
from that particular part of the orbit that today are still in the 
vicinity of the sun. In Fig.|2] we present the fraction of stars 
that can be found in the solar neighborhood as a function of the 
moment that the star escaped from the parent cluster 

The fraction of brothers and sisters of the sun that are still 
within lOOpc of our current location in the Galaxy is about 1% 
if the proto-star cluster dissolved shortly after the formation of 
the sun, and the fraction increases to ^ 8% if the cluster sur- 
vived longer (see Fig.|2]l. The fraction of siblings within a kpc 
is at least ^ 10 %. The majority is located around the past and 
future orbit of the sun in the Galactic disk and within about 
100 pc of the Galactic plane. This is a consequence of the rel- 
atively low velocity dispersion of the cluster compared to the 
orbital motion around the Galactic center. 

Today, we should still be able to recognize the siblings along 
the orbital trajectory of the solar system in the Galaxy. If the 
sun happens to be on a different orbit this will reflect in the 
orbits of its siblings. The distribution of proper motions and 
distances of these stars are presented in Fig. [3] where we show 
the result from a cluster of 2048 stars with a 1 pc virial radius 
that dissolved 4.6 Gyr ago. The proper motions and distances 
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of these stars change in a characteristic way, indicated by the 
soHd curve and in the direction of the two arrows. 

The majority of our lost sibHngs are easily identified by the 
Gaia astrometric satellite or by ground based searches. We are 
still surrounded by them, even though they hide among millions 
of other ordinary looking stars. What enables us to recognize 
the other siblings are their orbital characteristics, which should 
be comparable to the sun. In Fig.|4]we present the distribution 



of stars that once belonged to the proto-solar cluster; they are 
currently observable along the orbit of the sun in the Galaxy. 
The best place to look for the lost siblings is along the trajec- 
tory of the sun in the plane of the Galaxy, in leading and trailing 
orbits around the Galactic center. Identifying those stars will 
provide stringent limits on the sun's orbit around the Galactic 
center and gives us a unique window to study the conditions of 
the star cluster in which the sun was born. 
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Fig. 1 . — Constraints on the mass (M in Mq) and radius (R in pc) of tlie cluster in wliicli the sun was born. The allowed parameter space is fenced off by four 
curves, one for the minimum mass of about 500 Mq (dotted line), one limits the relaxation time (t^ix) to less than about 300 Myr (dashed curve), and two limiting 
the distance between which an encounter is likely to occur during the cluster lifetime (solid curves). If an enc ounter is expe cted outside renc = 10^ AU the orbits 
beyond Pluto are unlikely to have been affected, contrary to the observed high inclinations in the Kuiper belt (I da et alll2000!l . An encounter within rcnc = 100 AU 
is likely to be destructive to the solar system as we know it (Adams et al., 2006). In order to warrant one or more close supemovae near the sun before the cluster 
dissolves we require that the relaxation time is smaller than 300 Myr. The area fenced off by the dotted, the dashed and the solid curves is the favorable range of 
parameters for the proto- solar cluster. 
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Fig. 2. — The fraction of siblings of tlie sun which are still in the neighborhood as a function of how long ago the parental cluster dissolved. The stars within 
100 pc (bullets) and those within Ikpc (circles) are fitted by a straight line (dotted curves). The vertical error bars for the Ikpc calculations represent the one standard 
deviation from the mean and is based on calculating the trajectories of 2048 stars which were bom in a virialized Plummer sphere with M ~ 920 Mq and il = 1 pc 
in the Galactic potential. The fraction of nearby siblings depends only weakly on the uncertainty in the current position and velocity of the sun. 
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Fig. 3. — The characteristic pattern of the proper motion and distance of the stars in the solar neighborhood that once belonged to the same star cluster In this 
case we assumed that the star cluster dissolved quickly upon the formation and pollution of the sun, ~ 4.6 Gyr ago. The curve overplotted with the data points 
represent solar orbit with a slight offset of 1 pc along the x-axis and with a 1 Okm/s higher velocity in the x-direction. The two arrows indicate how the distance and 
proper motion change with time for stars that approach the solar position (right arrow) or are ahead of the sun (left arrow). 
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Fig. 4. — Top view of the Galaxy. The orbit of the sun for the last 4.6Gyr is presented with the thin dotted curve, starting at the arrow. The sun is currently 
located at a; = 8.5kpc with y = 0. A total of 1000 stars from a star cluster with the parameters used in Fig.|2]were followed together with the sun and evolved 
through time for 4.6 Gyr, their final positions are plotted as bullets. 



